The ultimate quantum limit to the linewidth of lasers 
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The standard quantum limit to the linewidth of a laser for which the gain medium can be 
adiabatically eliminated is £q = n/2n. Here k is the intensity damping rate and n the mean photon 
number. This contains equal contributions from the loss and gain processes, so that simple arguments 
which attribute the linewidth wholly to phase noise from spontaneous gain are wrong. I show that 
an unstimulated gain process actually introduces no phase noise, so that the ultimate quantum limit 
to the laser linewidth comes from the loss alone and is equal to l u \ t = n/4n. I investigate a number 
of physical gain mechanisms which attempt to achieve gain without phase noise: a linear atom- field 
coupling with finite interaction time; a nonlinear atom-field coupling; and adiabatic photon transfer 
using a counterintuitive pulse sequence. The first at best reaches the standard limit £q, the second 
reaches %to, while the third reaches the ultimate limit of £ u \ t = ilo. 



42.50.Ar, 42. 55. Ah, 42.50.Lc, 32.80.Qk 



I. INTRODUCTION 

It is more than 40 years since Schawlow and Townes in- 
troduced the idea of an "optical maser" ]lj] , now known of 
course as a laser. Probably the most famous result from 
this paper is the expression for the quantum-limited laser 
linewidth, their Eq. (17), 



2huj 

Pnut 



(Ac) 2 



(1.1) 



Here Aw osc is the half-width at half maximum (HWHM) 
of the laser, Alu is the HWHM of the relevant atomic 
transition, P out is the output power and lo the frequency 
of the laser. Defining £ = 2Aoj osc and 7 = 2Auj, this 
expression can be rewritten 
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huj 
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(1.2) 



where ST stands for Schawlow- Townes . The derivation 
of this expression assumes reabsorption of photons by 
atoms in the ground state of the relevant transition is 
negligible, and also ignores thermal photons and other 
extraneous noise sources. 

To accurately describe lasers a number of refinements 
must be made to the Schawlow- Townes expression [^). 
These are discussed in the Appendix. This discussion, I 
believe, helps to put in perspective some of the past work 
on quantum limits to the laser linewidth. The end result 
is that a better expression for the standard quantum limit 
to the laser linewidth is 



4t = 
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(1.3) 



Here st stands for standard (quantum limit). As ex- 
plained in the Appendix, ^baro is the bare linewidth, N 
is the number of coherent excitations stored in the laser 
mode and its gain medium and k is the cavity linewidth. 
The inequality is an equality only for perfectly efficient 
output coupling. 



In the limit 7 ^> k the gain medium can be adiabat- 
ically eliminated, resulting in Markovian evolution for 
the laser mode. This means that N can be replaced by 
n (the mean photon number), and £h&vc by k, to give the 
standard Markovian limit as 



to 
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(1.4) 



For the remainder of this paper I will assume the Marko- 
vian limit, and drop the adjective "Markovian" distin- 
guishing £q from £ st when no confusion is likely to arise. 

M ost older textbooks H-f| quote the result in 
Eq. ( |l.4| ), or one which reduces to it in the appropriate 
limit of neither reabsorption nor thermal photons. The 
first two of these HQ derive this result rigorously using 
quantum Langevin equations. All three attempt to ex- 
plain it in terms of the noise added by the spontaneous 
contribution to the (mostly stimulated) gain of photons 
from the atomic medium. Loudon S even recommends 
the argument based on the uncertainty principle given by 
Weichel §. 

The argument of Weichel is as follows (in my notation). 
In a laser at steady state, the rate of spontaneous emis- 
sions to total (spontaneous and stimulated) emissions is 
1 : n + 1. Since the total gain rate must equal the total 
loss rate *cn, the rate of spontaneous emissions is 



A = 



1 



(1.5) 



where it is assumed that n ^> 1. Now the reciprocal of 
this, At = 1/k, is || "the average time between phase 
fluctuations caused by spontaneous emissions into the 
mode." Invoking the uncertainty principle 



AEAt > h/2, 



(1.6) 



with the energy uncertainty of the mode being AE 
nhAuj osc gives 
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which agrees with the Schawlow-Townes result (1.2), 
with 7 replaced by k and P out by nn. 

The problem with all such simple arguments is that 
they put the blame for phase diffusion solely on the gain 
mechanism. This is an artifact of thinking in terms 
of normally-ordered operator products. That is, it re- 
sults from using (implicitly in most cases) the Glauber- 
Sudarshan P function [^[j9| as a true representation of the 
the fluctuations in the laser mode field. The P function 
is of course no more fundamental than the Q function || , 
which is a representation based on anti-normally ordered 
statistics. If one were to use the Q function as an aid to 
intuition, one would find that it is the loss process that 
is wholly responsible for the phase noise. Of course the 
rate of phase diffusion would agree with that from the 
P function, at least in steady-state where loss and gain 
balance. 

If one asks a question about phase diffusion, the only 
objective answer will come from using the phase basis it- 
self. This is far more difficult than using the more famil- 
iar phase-space representations, but some approximate 
results have been obtained pC| ]. These show that, at 
steady state, the phase diffusion has equal contributions 
from the loss and gain process. The same result occurs 
from a Wigner function calculation [Tc| ] . This is not sur- 
prising since symmetrically ordered moments are known 
to closely approximate the true moments for the phase 
operator for states with well-defined amplitude |ll[ . 

The fact that phase diffusion comes equally from the 
loss and gain processes suggests that the st and ard quan- 
tum limit to the laser linewidth, £q of Eq. (1.4), may not 
be the ultimate quantum limit. The contribution from 
the loss mechanism is unavoidable. A laser, at least in 
useful definitions | fL2| , requires linear damping of the laser 
mode in order to form an output beam. However, it may 
be that the standard gain mechanism could be replaced 
by some other gain mechanism that causes less phase dif- 
fusion. The ultimate quantum limit to the laser linewidth 
could thus be as small as one half of the standard limit. 

In this paper I investigate various gain mechanism in 
an attempt to find one which causes less phase diffusion 
than the standard gain mechanism. First, in Sec. II, 
I review the standard (ideal) model for a laser, giving 
rise to the standard quantum limit £q. Next, in Sec. Ill, 
I present gain without stimulated emission, which pro- 
duces a linewidth of £o/2, and discuss how this can be 
physically realized. In Sees. IV and V I present mod- 
els which attempt to approximate gain without stimu- 
lated emission, using a micromaser-like interaction and 
a nonlinear field-atom interaction respectively, and dis- 
cuss their success. After a comparison of these results in 
Sec. VI, I conclude in Sec. VII by returning to a deriva- 
tion of the quantum limits to the laser linewidth using 
an uncertainty relation. 



II. IDEAL STANDARD LASER 

A. The laser model 

The ideal standard laser master equation results from 
just about any gain medium under the appropriate (that 
is to say, ideal) conditions. Here I will present proba- 
bly the simplest derivation, in which the gain is due to 
the coupling of the laser mode with a single transition 
in an atom. Ignoring the other levels in the atom, the 
interaction Hamiltonian is 



H = in(aa^ - ah 



(2.1) 



where a is the annihilation operator for the cavity mode, 
a = \l)(u\ is the lowering operator for the atom, and f2 
is the one-photon Rabi frequency. 

Let the interaction time r be such that e\ffi <C 1, where 
e = fir and fi is the mean intracavity photon number. 
Then the unitary operator exp(— iHt) acting on the ini- 
tially factorized state R = p®\u){u\ can be expanded to 
second order in e to give the entangled state for the atom 
and field 

R = p <g> \u){u\ + e (a f p ® \l) (u\ + H.c.) 

+ e 2 (a f pa® |Z)(i| - \{aa\p\® \u){u\) . (2.2) 

Say there is a detector which detects which state the 
atom is in immediately after it has interacted with the 
field. If the outgoing atom is detected in the upper state, 
then the conditioned state of the field (the norm of which 
represents the probability of this detection result) is, to 
first order in e 2 , 



p u = (u\R\u) = (1 - e 2 A[a)\) p 

= exp(-e 2 aa t /2)/cexp(-e 2 aa' l '/2), 



(2.3) 



where the superoperator A[c] is defined for an arbitrary 
operator c by 



(2.4) 



If the atom is detected in the lower state (which happens 
rarely), the state is 

pi = (l\R\l) = e 2 J[a^}p, (2.5) 

where the superoperator J\c\ is defined by 

J[c)p = cpJ. (2.6) 

If this were all that there was to the model then the 
master equation would be found simply by averaging over 
the two results. Assuming that excited atoms enter the 
cavity as a Poisson process with rate r <C ft, the result 
would be 



p = Te 2 V[a)]p + nV[a]p. 



(2.7) 
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Here I have included linear loss (allowing the laser out- 
put) at rate k, and I am using the notation 



V[c] = J[c] - A[c\. 



As long as Te 2 < k, this master equation has a steady 
state. However it is not an appropriate steady state 
for the device to be considered a laser. As discussed in 
Rcf. jl2) , it is necessary to have n 3> 1 for the output of 
the device to be coherent (in a quantum statistical sense) . 
But i n th is limit, the stationary state of the master equa- 
tion (2.7) has a photon number uncertainty a{n) ~ n. 
This leads to enormous low- frequency (~ n/n) fluctua- 
tions in the intensity of the output beam. This ruins 
the second-order coherence of the device, which is why it 
could not be considered a laser 1131 . 



The origin of the problem with Eq. (2.7) is stimulated 
emission. Since this is part of the acronym l.a.s.e.r. it 
might be thought that this is a good thing, but it actually 
leads directly to t he enormous intensity fluctuations in- 
herent in Eq. ( |2.7| ). This is because stimulated emission 
implies that, for «>1, the intensity gain is proportional 
to the intensity. Thus if the intensity fluctuates above its 
mean value that fluctuation will be reinforced by an in- 
crease in the gain, and if it fluctuates below the mean 
the gain will correspondingly decrease. To avoid this, 
and hence obtain the output characteristic of a laser, one 
actually wants a photon gain which is independent of 
fluctuations in the intensity. 

In a real laser, this independence of gain rate with 
intensity occurs automatically as n becomes very large 
because of gain saturation. This is not difficult to derive 
in the master equation approach ]l3[ |. However the re- 
sult can be obtained more quickly in the context of the 
present model, by making sure that the atom gives up 
exactly one quantum of energy to the field, regardless of 
the field state. This is achieved by the following proce- 
dure. If the atom is detected in the lower state, then the 
field has gained a photon and the process can stop. If it is 
detected in the upper state, one must try again with the 
same atom (or, more realistically, another excited atom). 
This process continues until the atom is detected in the 
ground state. This process is shown in Fig. 1. 

Say K atoms are required before t he ( K + 1 ) th is de- 
tected in the lower state. From Eqs. ( ]2.3| ) and ( |2.5| ), the 
unnormalized state matrix after the (K + l)th atom is 

Pk = e 2 l /[a t ]exp[-ife 2 aa t /2]pexp[-ife 2 aa t /2]. (2.9) 

The norm of this state matrix is equal to the probabil- 
ity that this many atoms are needed. Thus, the average 
density operator, given that an atom is finally detected 
in the ground state, is 



(2.10) 



Using the fact that e 2 is small, the sum in Eq. ( 2.1G| ) 
can be converted to an integral by setting (3 — e 2 K: 



/>oo 

p' = J[(J]l cxp(-/3aa t /2) / 9exp(-/3aa t /2)d/3. (2.11) 
Jo 



(2 8) This can be formally evaluated |T^] as 



(2-12) 



The superoperator ^l[a^] _1 is well-defined because aa* is 
a strictly positive operator fl4|| . 

The action of the superoperator J^a^y^a^] -1 is to add 
a photon to the system irrespective of its initial state. 
That is to say, it shifts the photon number distribution 
upwards by one. If this addition of a photon is assumed to 
occur at Poisson-distributed times, with a rate T -C Hero, 
then a Markovian master equation for the field results. If 
one also includes linear damping at rate k as above, and 
lets the gain (the rate of photon addition) be T = up, 
then one gets 



kT 1 /? = p (jVMMP 1 - 1) P + T>[a]p 



= pV[a ] ]A[a 



tl-i 



V[a]p. 



(2.13) 



From Eq. (2.11) and the identity 



1= / du exp(— uaa^/2)aa ! exp(— uaa' /2) (2.14) 
Jo 

i t is easy to see that the ideal laser master equation 
( |2.13 ), first derived in Ref. is of the required Lind- 
blad form [ffll. 



B. Stationary State 



In the fock basis the laser master equation (2.13) is 



n + m 



Pn—l,m — l Pn,' 



n + m 



+ y/(n+l)(m + l)p n+ 



l,m+l- 



" Pn,m 

(2.15) 



Here, as in the remainder of the paper, I have set k = 1. 
Clearly the stationary state will be of the form p n>m = 
dn,m,Pn- The equation of motion for P n is 

K- x P n = p (P„_ x - P n ) + (n + l)P n+1 - nP n . (2.16) 

This has the stationary solution P n — e~^ l p n jn\. That is, 
the intracavity photon statistics are exactly Poissonian. 
The stationary state matrix can therefore be written 



\n)(n\. 



Equivalently, it can be written 

Pss = / d<t>\ae^){ae% 



(2.17) 



(2.18) 



where |a| = ^fp and \ae l< ^) is a coherent state of ampli- 
tude \ae 1 ^}. From either expression it is easy to verify 
that the mean number is Tr[o^ap S s] = p and the mean 
amplitude Tr[ap ss ] = 0. 



3 



C. Calculating the linewidth 



one obtains 



There are many different ways of calculating the 
linewidth of a laser from its master equation. One way is 
to covert the master equation into a Fokker-Planck equa- 
tion for a quasiprobability distribution function such as 
the P, Q or W function This is relat ively straight- 
forward for a master equation of the form (2.13), despite 
the apparent awkwardness of the inverse superoperator 
yl[at] _1 |l3). However, for other master equations as I 
will consider later in this paper, the conversion is not so 
simple. Therefore I will adopt a method using the Fock 
basis. The method is essentially a more rigorous version 
of that used by Sargent, Scully, and Lamb 

The linewidth I of a laser I have taken to be the FWHM 
of the Power spectrum 

f-OO 

P(u) oc / dr. 9 (1) (r)coswr, (2.19) 
Jo 

where the normalized first-order coherence function is 

3 «(r) = (ot(t + T)a(t)) ss /(a^a) ss . (2.20) 



If one represents the master equation ( [2.13 ) as p = Cp 
then one can write 



5 «(t) = Tr[af e^(ap ss )]fp 



(2.21) 



Note above that the stationary st ate m atrix p ss is a 
mixture of coherent states, as in Eq. ( 2.18 ). Since g^(r) 
is invar iant u nder a phase shift, Eq. (2.18) implies that 



in Eq. ( 2.21 ) one can take p s 
Then Eq. fl2.21|) becomes 



| a) (a |, with \a\ 



gW(r) = Tr[J ap(r)}/p, 
where p(t) obeys the master equation ( |2.13 ) and 
p(Q) = \a)(a\. 

If one defines 

/„(*) = y/np n -i tn (t)/a* 
then one can write 

5 (1) w = E/«w- 



(2.23) 



(2.24) 



(2.25) 



Clearly if one can determine the evolution of f n {t), one 
can f ind g ^> (t) and hence the linewidth of the laser. From 
Eq. (|2.13D one finds 



fn = A* 
Defining 



2)i 



2n-l 



fn-l - Hfn + nf n+ i 



2n- 1 



r„(t) = 



n.fn+i{t) ' 



fn- (2-26) 



(2.27) 



fl I 



2n(n — 1) fi 2n — 1 



/„■ (2.28) 



Now from the definition ( |2.27 ), r n (0) = 1. Assum- 
ing that this ratio remains unity, expand Eq. ( 2.2S ) to 
leading order in 1/p, to get 



fn 



An 



fn ■ 



(2.29) 



Solving this and substituting into Eq. (2.27) gives, to 
leading order, 



t 



t 



(2.30) 



where the expansion to first order is valid for times much 
less than p 2 . Since, as will be shown, the coherence time 
~ 2/t is of order p, it is quite safe to make this expansion 
even for times long compared to the coherence time. 

Substituting this expression for r n (t) into Eq. ( 2.28| ) 
gives the more accurate expression 



1 

4n 



Since the initial condition is 



fn • 



MO) = e-» 



/■'" 



(n-1)!' 



(2.31) 



(2.32) 



(2.22) th e only significant contribution to the sum ( [2.25] ) comes 
from n such that |n— p\ < JJi. Also, as noted abov e, one 



can assume t < n. Then the correction term in Eq. (2.31) 
is of order /i -1 / 2 and can be ignored. One can thus re- 
turn to the expression Eq. (2.29), which becomes (again 
ignoring corrections of order p / 2 ), 



fn 



1 

Tp 



fn 



The first order coherence function is thus 
5 «(t) = exp(-r/4/x), 



(2.33) 



(2.34) 



so that the coherence time is 4/x (which is of order p as 
promised). The Fourier transform of this expression is a 
Lorentzian with FWHM 



1 

2p' 



(2.35) 



This is the standard quantum limit Iq of the linewidth 
for an ideal laser. 
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III. GAIN WITHOUT STIMULATED EMISSION 

Since the "stimulated emission of radiation" is part of 
the acronym for laser, it might be thought that stimu- 
lated emission is essential to produce a laser. While a 
typical laser does rely upon stimulated emission to en- 
sure th at it runs single-mode, the fact that the model of 
Section II A adds photons one by one suggests that it is 



not strictly necessary. I will now show that stimulated 
emission is indeed not necessary for laser action, and in 
fact that eliminating stimulated emission eliminates the 
phase diffusion caused by the gain process. 

Stimulated emission is a simple consequence of th e lin - 
ear coupling of the laser field to it s so urce, as in Eq. (2.1). 
That is to say, the Hamiltonian (2.1) is linear in the an- 
nihilation operator a which, for classical fields, can be 
replaced by the c- number y/ne 1 ^. Whenever a rate is cal- 
culated in quantum theory it depends on the square of 
the Hamiltonian. Hence the fundamental gain rate from 
a linear coupling will vary as n, which is the so-called 
stimulated emission or Bose-enhancement factor. A fully 
quantum calculation of course gives spontaneous emis- 
sion as well, and hence a gain rate proportional to n + 1. 

Since stimulated emission can be traced to the pres- 
ence of a in the coupling Hamiltonian, the only way to 
totally remove it is to substitute for a a different lowering 
operator, one whose classical a nalo gue does not increase 
with n. That is to say, in Eq. (2.1), replace 



= ^2 Vn\n- l)(n\ 



(3.1) 



by the Susskind-Glogower [hjj e = e 1 ^ operator 



(aa 



l)(n| 



(3.2) 



71=1 



The new Hamiltonian would be extremely nonlinear if 
expressed as a power series in a and , but it cannot be 
denied that it will not exhibit any stim ulate d emission. 

Replacing a by e in the Hamiltonian (2.1) presents no 
problems in the rest of the derivation. Moreover, it is 
not even necessary to assume that e = Qt is very small. 
Instead, the result is independent of e, due to the fact 
that ee^ = 1. In particular, if one chooses e = 7r/2, the 
transformation effected on the field by one transit of the 
atom is semi-unitary: 



exp 



Here \tb) is the state of the field and 



S = et = J2\n + l)( 



(3.3) 



(3.4) 



n=0 



The operator S is semi-unitary rather than unitary be- 
cause St£ = 1 but SS* = 1- 



Surprisingly, this transformation can be achieved phys- 
ically using only the usual electric-dipole coupling p8| . 
The trick is to use a three-level A atom and another, clas- 
sical field ps|| . Then, using a using a counter- intuitive 
pulse sequence, the atom is transferred from one lower 
state to the other, and one photon is created in the cav- 
ity field (with the energy lost from the classical field). 
Like the gain process in Sec. II, this adds precisely one 
photon to the field. The difference is that it does this 
without entangling the state of the field and the atom, 
and hence leaves the state of the field pure. 

Taking the rate of addition of ph otons to the field to 
be r as before, in place of Eq. (2.13) one obtains 



p = nV[eJ}p + V[a]p. 



(3.5) 



In the fock basis this becomes 



= P(p n -l,m-l - Pm,n) ~ (n + m)p, hm /2 



y/(n + l)(m+ l)p„+l,m+l- 



(3.6) 



This yields exactly the same equation for the diagonal 
elements (the photon number populations). Hence the 
unstimulated master equation produces exactly the same 
photon numb er sta tistics as does the standard laser mas- 
ter equation ( 2.13 ). 

To calculate the linewidth, proceed as before. One 
finds the following equation for f n , defined as in Sec. II: 



fn = V 



fn-1 - fn) + nf n+1 - 2n - f n (3.7) 



n- 1 



\Jn(n - 1) r„_i — p -\ 

r. 



p In — 1 



(3.8) 



Assuming r n « 1 yields, as above, the self-consistent so- 
lution 



fn 1=3 fn- 

8p 

The first order coherence function is therefore 
5 W(r) = exp(-T/8/*), 



so that the linewidth is 



1 



(3.9) 



(3.10) 



(3.11) 



This is half the standard quantum limit £o of Eq. (1.4 ). 
As explained in the introduction, the standard quantum 
limit for the phase diffusion rate contains equal contribu- 
tions from the gain and loss processes. The gain process 
considered in this section does not introduce any phase 
noise; the operator is more or less the exponentiation 
of the phase operator and so increases the photon num- 
ber without affecting the phase distribution at all. Thus 
the phase diffusion in this model comes wholly from the 
loss process, and the rate is half the standard rate. 
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IV. FINITE ATOM-FIELD INTERACTION TIME 

The preceding section showed that an interaction in 
which the atom is sure to give up its quantum of energy 
to the field from a single pass results in a linewidth a fac- 
tor of two smaller than the standard limit. It was noted 
there that this could be achieved using an adiabatic pas- 
sage, but this has yet to be done experimentally. This 
suggests that it would be worth exploring other ways to 
mimic the unstimulated gain process. 

In this section I investigate one idea, based upon the 
gain mechanism of a micromaser [20 2M. This utilizes 
the same Jaynes-Cummings coupling ( |27l| ) as in Sec. II. 
The difference is that the scaled interaction time e = Qt 
is not assumed to be small. This modifies the results of 
Sec. II as follows. The state of the field conditioned on 
the detection of an atom in the lower state is lEOI] 



Pi = Jip, 



where 



Ji = J 



e' sin 



(4.1) 



(4.2) 



The field state conditioned on an atom passing through 
and remaining in the upper state is |20l 



Pu = JuP, 



where 



cos 



(4.4) 



For states having a photon distribution localized 
around n, if e is such that e^/fi w n/2, then it would 
seem that the action of the above superoperators could 
be approximated by 



Ji 

Ju 



0. 



(4.5) 
(4.6) 



That is, the atom would almost certainly come out in the 
lower state, having given up its quantum of energy to the 
field. This is the same situation as for the unstimulated 
gain as shown in Sec. III. This is why a finite interaction 
time e might be expected to lead to a linewidth below 
the standard limit. 

If atoms are injected at a Poissonian rate /i then the to- 
tal master equation is the usual micromaser master equa- 
tion 



p = {p{J u + Ji-l) + -D[a]}p. 



(4.7) 



Here linear damping at rate unity has been included also. 
This master equation has very complicated dynamics. 
For some values of e and ^ the stationary state does not 
have a well-defined intensity. That is, it is not the case 
that u(n) <C n. Hence the device is not necessarily a true 
laser in the sense of Rcf . wA . 



To ensure that the a well-defined photon number dis- 
tribution is produced, the same technique as in Sec. II 
can be used. That is, if an atom is detected still in the 
upper state it is sent through again until it is detected in 
the lower state. The resulting master equation is 



\nJiY,JZ + V[a]\p 

I k=0 ) 

[nJiil-J^+Vla]) 



(4.8) 
(4.9) 



In the photon number basis 

sin(e v / n) sin(ey / ro) 



Pn,: 



1 — cos(e^/n) cos(ey/m) 

- PPn.rn ~ (n + m)/9 n , m /2 

+ v/(n+l)(m + l) 



(4.10) 



To find the linewidth, one proceeds as before to get the 
following equation for f n : 



fn = P 



3in(evn— T) sm(e^/n) 



y/n — 1 [l — cos(eV« — 1) cos(ey / n)] 

■ Pin + Tifn+1 Z Jn- 



fn-1 



(4.3) Using the parameter 



(4.11) 



(4.12) 



one can continue the analysis as before and find eventu- 
ally 



1 



/x 2 sm 2 ((j)/n) 
sin 2 d> 



fn 



(4.13) 



That is, the linewidth of the laser is found to be 

2" 



1 

~P 



1 



sm(<f>/p) V 
(sin <t>)/n) 



(4.14) 

It is easy to verify that this expression has a global 



minimum 



lim 



1 



1 



sin(<ft//z) V 
(sin <f>)/[i) 



1 

2/? 



(4.15) 



The limit ip — > is the limit of short interaction times in 
which the original model of Sec. II is recovered, and the 
original linewidth £q also. That is, no linewidth narrow- 
ing is possible using a finite interaction time in preference 
to an infinitesimal interaction time, despite the fact that 
the former can deposit a photon in the cavity in a single 
pass of the atom with very high probability. 

This line-broadening is definitely not an artifact of the 
assumption that the atom is always put through again if 
it is detected still in its upper state; a similar result is 
obtained for the usual master micromaser equation with 



G 



a single pass per atom p2fl. The approach to calculating 
the linewidth used in Ref. p2| was similar to the one used 
here. A more accurate estimation of the linewidth for the 
usual micromaser has to take into account the fact that 
the intensity is not always well defined ^3|. This yields 
some deviations from the simple theory of Ref. |2^j , but 
still never shows any line-narrowing. 

The reason that no linewidth narrowing occurs can be 
seen from the method of calculation I have employed. 
What turns out to be crucial is not to try to mimic the 
two terms in the ideal unstimulated gain, namely 



V[e*]=J[et]-A[J]=J[J]-l, 



(4.16) 



but rather to mimic the following ratio of matrix elements 
involving these two terms: 



(n-l\{j[J]\n)(n + l\}\n) 
(n|{^[et]|n)(n+l|}|n + l) 



1. 



(4.17) 



In the unstimulated case the ratio is unity, and the dif- 
ference from unity in other cases is proportional to the 
contribution to the linewidth from the gain process. For 
the standard ideal laser, 



(rc-l|{,7[at]|n)(n + l|}|n) 
(n|{>t[at]|n)<n + l|}|n + l) 



(4 - 18) 



Multiplying the deviation from unity by the gain constant 
fx and replacing n by the mean photon number \x gives 
1/8/j.. This is the standard contribution to the linewidth 
from the gain. For the micromaser, 



(rc-l|{Ji|n)(n + l|}|n) 
(n|{[l- J u ]\n){n + l\}\n + l) 



1 



sin 2 (^)//i) 



(4.19) 



which again explains the result in Eq. ( 4.14 ). 



V. NONLINEAR ATOM-FIELD INTERACTION 

With now a better understanding of how to reduce the 
gain-induced phase diffusion, I turn to a second method 
for trying to mimic unstimulated gain. As noted in 
Sec. Ill, the operator e* would require an infinite series 
to be expressed in terms of powers of a and aJ. Any 
Hamiltonian containing infinite powers of the field is un- 
physical. However, nonlinear optical processes containing 
field powers greater than unity do occur. This suggests 
that it is worth considering the following approximation 



* = at(aat)- 1 / 2 = a) [/, + {aa} - »)] ~ 1/2 
at /3 laat N 
^ \2~2~ix~ 



(5.1) 
(5.2) 



That is, I wish to consider a nonlinear Jaynes-Cummings 
Hamiltonian of the form 



(3 - aa f In) + (3 - aa) / p)aa^] , (5-3) 

which I expect to be useful when the photon number is 
approximately p. 

Physically, this Hamiltonian means that there are two 
processes which can excite the atom. The first is the 
usual linear dipole coupling to the field. The second is 
a three-photon process whereby a photon is virtually ab- 
sorbed and re-emitted before finally being absorbed by 
the atom. The Hamiltonian matrix element for the sec- 
ond process is much smaller (for p ^> 1), which is phys- 
ically reasonable, and is of the opposite sign. It is very 
doubtful that such a Hamiltonian could be achieved sim- 
ply using a two-level atom. However, it is possible that 
an effective Hamiltonian of this form could be achieved 
using a multilevel atom, and other fields. I will not fur- 
ther discuss the feasibility of producing this Hamiltonian, 
as my chief concerned is with the questi on o f principle: 
how well can the nonlinear Hamiltonian (5J3) reproduce 
the results of the unstimulated laser? 

Assuming, as in previous sections, that the atoms are 
initially in the upper state and that any atom which exits 
the cavity still in the upper state is put through again, 
one derives, following the method of Sec. II, the following 
master equation for the cavity mode. 

p = pV[a){Z - fla* / 'p)]A[a) (3 - aa,y fi)}- 1 p + V[a}p. 

(5.4) 

This has the same Poissonian mixture of number states 
as in the standard laser, and is amenable to the same 
method of calculating the linewidth. The result is 



3 

8/? 



(5.5) 



That is, the contribution from the gain is 1/8/x, which 
is half the standard result and half the contribution of 
l/4/i from the loss (which is of course unchanged). This 
result can again be understood from the ratio 



(n - 1| {JV(3 - aaVp)}\n) (n + 1|} \n) 
(nj MM(3 - aa^/p)}\n)(n + 1|} \n + 1) 



1 - 



1 



16n 2 
(5.6) 



VI. DISCUSSION 

The standard quantum limit to the laser linewidth is 
not the ultimate quantum limit, even for the Markovian 
case in which the gain medium is eliminated from the 
equations of motion of the laser mode. Hidden within 
the standard Markovian expression 



(6.1) 
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are equal contributions of k/Au from the gain and loss 
mechanisms for the laser. The latter contribution is a 
fundamental limit because linear loss is necessary for a 
coherent output beam to form. However the former re- 
sults from a particular (extremely reasonable) assump- 
tion about the gain mechanism for laser action, namely 
that it comes from a weak linear coupling between the 
field and the gain medium. 

These arguments suggest that a different sort of gain 
mechanism could produce a laser with a linewidth up to 
50% below the standard quantum limit. As I have shown 
above, this ultimate Markovian limit 



4it - 



Art 



(6.2) 



can be achieved with a gain mechanism in which stim- 
ulated emission into the laser mode is eliminated. This 
requires that the matrix element for the addition of a 
photon to the laser mode be independent of the number 
of photons in the mode. As discussed, this could be phys- 
ically achieved with adiabatic transfer of photons from 
another field using a counter-intuitive pulse sequence. 

I also examined two other gain mechanisms with sim- 
ilarities to the non-stimulated gain, to see if they also 
produced linewidth narrowing. The first, using the usual 
Jaynes-Cummings Hamiltonian but with a finite interac- 
tion time (as in the micromascr) did not. The second, us- 
ing a nonlinear Jaynes-Cummings Hamiltonian involving 
three-photon as well as one-photon processes, produced 
a linewidth of 



t = 



3k 



(6.3) 



That is, the phase diffusion due to the gain was reduced 
by 50% from the standard limit, resulting in an overall 
reduction of 25% in the linewidth. Presumably higher- 
order nonlinear optical processes could more closely ap- 
proach the ultimate limit. However, the difficulty in pro- 
ducing such nonlinear optical processes, and the fact that 
even a third-order nonlinearity fails to reach the ultimate 
limit, suggests that the adiabatic transfer is a better ex- 
perimental option for probing the ultimate quantum limit 
to the laser linewidth. 

The ultimate limit for the rate of phase diffusion at- 
tained by eliminating gain noise can also be obtained, for 
short times, by instead eliminating loss noise. This can 
be achieved by coupling the laser output into a squeezed 
vacuum rather than a normal vacuum |^,^]. This only 
works for short times because it requires a specific phase 
relation between the squeezed vacuum and the coher- 
ent field in the laser, which will not remain valid since 
the laser phase continues to diffuse. It was suggested in 
Ref. [g3J that it might be possible to produce the squeezed 
vacuum by driving the squeezing device with the laser it- 
self. In this case the whole squeezing device should really 
be considered as part (an internal absorber, in fact) of 
the laser, so that n in the original laser cavity should 
no longer be used as a good measure of the total stored 



excitation. Similar comments could be made about the 
proposal of Ghosh and Agarwal |2(|, who also misquote 
the expression for the standard quantum limit given in 
Ref. M by a factor of two as their Eq. (18). I believe 
that a rigorous analysis of these proposals would reveal 
no reduction below the standard quantum limit. 



VII. CONCLUSION 

In the introduction I reproduced a simple argument 
purporting to use the time-energy uncertainty principle 
to derive the standard laser linewidth as a consequence 
of phase diffusion due to the gain process. The results 
of this paper show that any such simple argument is un- 
tenable since the gain process contributes only half of 
the standard phase diffusion rate. To compensate for 
disposing of this simple argument, I will conclude this 
paper with a (not quite so simple) argument deriving the 
ultimate Markovian quantum limit £ u \t from another un- 
certainty principle argument. 

Instead of the energy-time relation, which is of doubt- 
ful content, I will use the quadrature uncertainty relation 



V(X)V(Y) > 1, 



(7.1) 



where X/2 and Y/2 are the real and imaginary compo- 
nents of the laser mode amplitude a. Clearly the vacuum 
state is rotationally symmetric with 



V(X) = V(Y) = 1, 



(7.2) 



and this holds also for a coherent state (which is the state 
the laser mode can be assumed to be in). 

Let the mean amplitude of the coherent state be real 
and positive so that X/2 = y/E and Y = Q. The phase 
variance is 



V{4>) = V ( artan— 



von 

X 2 



for ii » 1. Now the effect of linear damping for an in- 
finitesimal time dt is to reduce the mean photon number 
of the coherent state from n to n(l — Kdt). Thus the 
change in the phase variance is 



(7.4) 



A noiseless gain process will return the mean photon 
number to n without increasing the phase noise. There- 
fore the phase variance increases at least as 



m ~ | . 



(7.5) 



The linewidth is defined from the two-time correlation 
function 

(a\t)a(0)) ~ n(e^ (t) ) ~ ne~ vw/2 ~ ne~ Kt/m . (7.6) 
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The Fourier transform of this expression is a Lorentzian 
with a FWHM of 



K 



(7.7) 



which is the ultimate quantum limit to the laser 
linewidth, as claimed. 
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APPENDIX: REFINING THE 
SCAWLOW-TOWNES LIMIT. 



The Schawlow-Townes expression 



-7 



(Al) 



was derived in the days before good optical cavities, and 
hence implicitly assumes that the atomic linewidth 7 
is much smaller than the (FWHM) cavity linewidth k. 
With k < 7, it is necessary to replace 7 by the bare 
linewidth of the laser ibarc- This is the frequency spread 
the output would have if the pump were suddenly turned 
off and all of the energy allowed to escape. For a large 
class of lineshapes, it can be shown that a reasonable 
approximation to the bare linewidth including contribu- 
tions from the atomic (or other gain) medium and the 
cavity is 



^barc I 



(A2) 



For instance, this expression agrees with that given by 
Haken (p. 103 of Ref. citeHak84) for the case where 
k > 7. In the other cases, where k <C 7, the 7 in 
the Schawlow-Townes expression is sim ply replaced by 
k which also agrees with Eq. (|A^). The corrected 
Schawlow-Townes expression is thus 



e' ST — 



9 l 

' bare 



7 2 k 2 



(7 



(A3) 



The second correction which must be made to the 
Schawlow-Townes linewidth relates to its use of the out- 
put power. Say, for argument's sake, that one one has 
a laser with linewidth given by the Schawlow-Townes 
limit, with all of the power coming out of one mir- 
ror. Then say that the mirror is replaced by one of 
the same reflectance, but with larger internal absorp- 
tion. Then the power loss per round trip is identical, 
so the laser dynamics remain the same and the linewidth 
would remain the same. But the power out would be re- 
duced because the transmittancc is reduced. Therefore 



the Schawlow-Townes formula would now predict an in- 
creased linewidth, which does not occur. In other words, 
the actual new linewidth would be less than the quantum 
limit set by the Schawlow-Townes formula. It is obviously 
inappropriate that a quantum limit can be surpassed by 
building a worse device. 

The resolution to this problem with the Schawlow- 
Townes linewidth is to eliminate P ou t from the expression 
by recognizing that 

(A4) 



^barc 

is an upper bound on the mean energy E stored as co- 
herent excitations in the laser system. If all of the stored 
coherent excitation eventually makes it into the output 
beam of the laser then the bare linewidth £bare is due 
wholly to the output coupling and P out = ^bareP- In gen- 
eral Pout is less than this. Reducing the output coupling 
efficiency (as discussed in the preceding paragraph) will 
not affect E so it is the correct parameter to use, rather 
than P out . The doubly corrected Schawlow-Townes limit 
is thus 



„// _ ^barc^W 
^ST — Jf, S 



flLU 



7 2 K 2 



Pout (7 + «) 2 



(A5) 



where the inequality becomes an equality only for per- 
fectly efficient output coupling. 

It is convenient to define the number of quanta of co- 
herent excitation, N — E/T110. For the case k ^> 7 the 
excitation stored in the gain medium is negligible and 
N = n, where the latter represents the mean photon 
number in the cavity. If the gain medium cannot be 
adiabatically eliminated then N must include the exci- 
tations stored coherently in the gain medium as well. If 
7 C k, as in the original Schawlow-Townes expression, 
these excitations in the gain medium will be the domi- 
nant ones. 

The final correction which needs to be made to the 
Schawlow-Townes linewidth is to insert a factor of | . The 
Schawlow-Townes limit without this factor is appropriate 
to a laser below threshold in which the complex ampli- 
tude of the field undergoes large slow fluctuations (for 
N 3> 1 which is the limit in which the Schawlow-Townes 
equation is valid). Above threshold, the laser intensity 
fluctuations are almost eliminated , leaving only phase 
fluctuations. This increases the coherence time by a fac- 
tor of two, so that the final corrected expression for the 
laser linewidth is 



2N 



< 



7 2 K 2 



2Pout (7 + *0 2 ' 



(A6) 



Here st stands for standard (quantum limit) as opposed 
to ST which stands for Schawlow-Townes . 

In the limit k <C 7, which is the usual limit for most 
lasers, and which allows the gain medium to be adiabat- 
ically eliminated from the field equations, one obtains 
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K ^ TliO 

0= 2fi- 2P out ' 



(A7) 



This result has often (including by myself [Q) been 
quoted as the Schawlow-Townes limit, despite the ob- 
vious differences from Eq. (Al). Here I will call it in- 



stead the standard Markovian quantum limit to the laser 
linewidth. "Markovian" refers to the fact that the equa- 
tions of motion for the laser mode, including gain and 
loss, are well-approximated by Markovian equations. For 
the gain process this is a consequence of adiabatically 
eliminating the gain medium. For the loss process, it is 
simply a consequence of assuming a high-Q cavity. Cor- 
rections (upwards) for non-Markovian loss (low-Q cavi- 
ties) are discussed for example in Ref . (^?J , but here I will 
always assume a high-Q cavity. 

Obviously for 7 < n, the linewidth of Eq. ( [A6| ) will 
be le ss t han the standard Markovian quantum limit of 
Eq. (A7). That is a reflection of the fact that in this 
case the bare linewidth £bare is less than k, and also that 
the gain medium is an extra reservoir of energy (coher- 
ent with the laser mode) so that N is greater than n. A 
linewidth which, in the ideal limit, reduces to 4t was re- 
cently derived in Ref. 2g|], for a laser with 7 < k. These 
authors claimed that this was "reduced compared to the 
Schawlow-Townes limit" because they followed the com- 
mon (but, in my opinion, erroneous) practice of iden- 
tifying Iq as the Schawlow-Townes limit. To me this 
seems to be an example of imprecise terminology obscur- 
ing an otherwise valuable contribution to fundamental 
laser physics. In this paper I always work with models in 
which the gain medium can be adiabatically eliminated, 
so that £ = 4t- 
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State 
Detector 



FIG. 1. Schematic of an ideal laser. An atom in the upper 
state passes through the cavity and its state is then detected. 
If the atom remains in the upper state, the process is repeated 
until it is detected in the lower state. The time for this process 
(including repetitions) is assumed to be very short compared 
to the cavity damping time kT 1 . Once the atom is detected 
in the lower state, a new upper-state atom is injected after 
a random waiting time r having an exponential distribution 
w(t) — cxp( — k,(it). Here /i is the desired mean number of 
photons in the cavity. 
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